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Abstract —  This  paper  expands  on  the  author’s  previous  work 
by  adapting  the  CLEAN  algorithm  to  address  low  signal  to  noise 
(SNR)  targets.  The  Reformulated  CLEAN  Detector  is  presented 
which  is  shown  to  allow  the  detection  of  low  SNR  targets  in  the 
presence  of  large  targets.  Performance  results  are  presented 
that  show  how  good  performance  can  be  attained  by  combining 
the  correlator,  CLEAN  Deconvolver  and  Reformulated  CLEAN 
Detector 

I.  Introduction 

The  CLEAN  algorithm  was  developed  by  astronomers  [1] 
to  enable  them  to  image  dim  objects  in  the  presence  of  bright 
objects.  The  problem  that  CLEAN  tries  to  overcome  is  that 
the  bright  or  high  amplitude  objects  have  responses  that 
spread  out  and  obscure  dimmer  objects.  CLEAN  assumes  that 
the  response  to  a  point  object  is  known  and  is  called  the  point 
spread  function  (PSF).  It  is  further  assumed  that  the  largest 
amplitude  response  in  the  scene  is  the  location  and  amplitude 
of  a  valid  object  (as  opposed  to  a  spurious  response  due  to 
noise  or  the  sidelobes  of  other  objects).  The  largest  amplitude 
and  its  location  is  used  to  position  and  scale  the  PSF  such  that 
the  largest  object  can  be  subtracted  and  thus  removed  from  the 
scene.  This  will  reveal  previously  obscured  objects.  This 
process  is  repeated  until  all  that  is  left  in  the  scene  is  noise 
(and  the  unavoidable  errors  in  the  CLEAN  algorithm).  The 
previously  noted  peak  amplitudes  and  their  positions  are  taken 
as  the  valid  objects  in  the  scene.  Thus  the  CLEAN  algorithm 
is  an  approach  to  deconvolution  where  the  PSF  is  removed 
from  the  image.  There  are  a  number  of  papers  applying 
CLEAN  to  astronomy  with  [2]  providing  a  recent  overview  of 
the  application  of  CLEAN  to  astronomy.  The  CLEAN 
algorithm  is  also  of  interest  to  the  radar  community. 
Consequently,  there  are  a  number  of  papers  proposing 
modification  to  CLEAN  and  its  application  to  radar  problems 
[3-8]. 

This  paper  extends  the  results  of  the  author's  two  previous 
papers  on  application  of  the  CLEAN  Algorithm  to  the 
condition  of  low  signal  to  noise  ratio  (SNR)  targets.  The  first 
paper  [9]  reviewed  the  CLEAN  algorithm,  attempted  to  put  it 


on  a  more  theoretical  ground  and  derived  the  CLEAN 
Deconvolver  and  the  CLEAN  Detector.  The  CLEAN 
Detector  was  derived  using  the  Neyman-Pearson  Theorem 
and  provides  the  maximum  probability  of  detection  for  a 
specified  probability  of  false  alarm.  The  ambiguity  function 
showed  extremely  good  performance  and  the  mean  square 
error  of  estimating  the  target  scatterers’  amplitudes  were 
shown  to  be  significantly  better  than  the  correlator.  The 
second  paper  [10]  addressed  the  performance  of  these 
processors  for  the  condition  of  Doppler  mismatch.  It  adapted 
the  CLEAN  Deconvolver  and  the  CLEAN  Detector  to 
unknown  Doppler. 

For  the  current  work  the  CLEAN  Deconvolver  will  be 
taken  as  the  starting  point.  The  CLEAN  Deconvolver  gets  its 
name  from  the  fact  it  is  the  straight  forward  deconvolution  of 
the  received  data  to  recover  the  scene.  This  can  be  seen  with 
the  baseband  signal  model  as 


y  =  WHc  +  n  (l) 

where  multiplication  by  W  performs  the  convolution 
operation  and  is  structured  as 


0 

9 

0 

vector  of  the  baseband 
version  of  the  transmitted  signal,  c  is  the  vector  of  the  target 
complex  in  reverse  range  order  and  n  is  Gaussian  noise.  Note 
that  it  is  presumed  that  there  is  a  scatterer  in  each  range  cell 
with  unknown  complex  amplitude. 

The  CLEAN  deconvolver  given  in  [9]  is 

c  =  (WWfl)'1  Wy.  (2) 


0 
a' 

|_wf  •••  0 

w  is  the  complex  conjugate  of  the 
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Some  authors  call  this  the  least  squares  estimator  [11]  which 
it  is.  However,  it  is  important  to  note  that  this  is  an  estimator 
that  conforms  to  the  linear  model  and  is  therefore  based  on 
theorem  4.1  of  [12,  page  85-86],  the  Minimum  Variance 
Unbiased  (MVU)  Estimator  for  c.  This  means,  for  example, 
that  (2)  will  give  a  smaller  mean  square  error  in  estimating  c 
than  the  matched  fdter  even  though  the  matched  fdter  gives 
the  maximum  SNR  for  same  problem.  (This  is  verified  in 
[9].)  The  advantage  of  using  (2)  as  a  starting  point  for 
CLEAN  is  that  its  performance  is  optimal,  it  is  easy  to 
calculate  and  gives  a  lower  bound  for  the  estimation  error. 

II.  Signal  Model 

The  problem  that  this  paper  addresses  is  detecting  small 
targets  in  the  presence  of  a  large  target  under  the  condition  of 
low  SNR.  In  order  to  remove  much  of  the  artificiality  of  a 
baseband  only  analysis  this  paper  uses  a  simulation  of  a 
modulated  carrier  pulse.  This  gives  a  signal  model  as 

x(t)=  s{t)cos{2rfc)*Yjai5{t-ti)  +  n{t)  (3) 

1=1 


where  *  is  convolution,  s(t )  is  the  baseband  pulse,  fc  is  the 
carrier  frequency,  q  is  the  number  or  targets  (or  scatterers),  a, 
is  the  complex  amplitude  of  the  target,  6(  )  is  the  Dirac  delta 
function,  t,  is  the  range  induced  time  delay  and  n(t)  is  white 
Gaussian  noise.  Figure  1  shows  the  standard  receiver 
processor  that  produces  the  data  vector  y  based  on  the 
returned  signal  x(t). 


-2sin(2n fct) 

Equivalent  Processing  Block  Diagram 

Figure  1. 


In  order  to  do  a  digital  simulation  of  Figure  1,  a  Matlab 
simulation  was  coded  that  implemented  the  processing  in 
Figure  2.  In  this  simulation  the  initial  sample  rate  was  set  at 
five  times  (M= 5)  the  ultimate  sample  rate  required  for  the 
Nyquist  sampling  of  s(t).  This  allowed  for  setting  target 
locations  at  factional  sample  intervals. 


Ts  7 ;  =T'tM 


2exp(-\2nfckTs) 


Matlab  Simulation  Block  Diagram 

Figure  2 

Now  it  is  instructive  to  relate  Figure  2  and  (3)  to  (1)  and 
(2).  The  PSF  is  w  and  is  equal  to  the  output  y  of  Figure  2 
when  the  input  x(t)  is  simply  s(t)cos(2nfcf),  n  is  equal  to  the 
output  y  when  x(t)  =  n(t),  and  the  elements  of  c  are  ck  =  al; 
where  k  =t/Ts  rounded  to  the  nearest  integer. 

III.  Reformulated  CLEAN  Detector 

While  the  CLEAN  Deconvolver  is  the  preferred  processor 
to  detect  and  estimate  closely  spaced  targets,  it  can  have 
significant  performance  problems  in  situations  of  low  SNR. 
In  this  case  the  preferred  processor  is  the  matched  filter, 
however  the  range-time  sidelobes  can  interfere  with  the 
detection  of  smaller  targets.  The  approach  to  overcome  these 
problems  is  to  develop  a  new  detector  structure  that  combines 
the  benefits  of  the  CLEAN  Deconvolver  and  the  matched 
filter. 

The  original  CLEAN  Detector  was  shown  to  be  equivalent 
to  a  deconvolver  [9].  It  also  handled  the  situation  of  eclipsed 
targets.  This  new  version  of  the  CLEAN  Detector  will 
assume  that  there  are  no  eclipsed  targets  (as  the  CLEAN 
Deconvolver)  and  will  simultaneously  act  as  a  matched  filter 
for  small  targets  while  eliminating  large  targets. 

To  begin  with  it  is  assumed  that  all  large  targets  (i.e., 
targets  large  enough  to  suppress  small  targets  with  their 
sidelobes)  have  been  detected.  This  can  be  done  either  by 
using  the  CLEAN  Deconvolver  or  the  matched  filter.  The 
amplitude  and  position  of  these  targets  are  used  to  estimate 

C+  which  will  contain  the  amplitude  of  large  targets  and  zero 
for  all  other  elements.  To  derive  this  new  detector  the  signal 
part  of  y  for  a  target  in  the  Cth  range  cell  model  is  defined  as 

ys=W8kak  (4) 

where  8k  is  all  zeros  except  for  the  A'th  element.  The 
interference  part  of  y  (large  targets  and  noise)  is  given  by 

y;=WV+n  (5) 

Based  on  (5)  the  correlation  matrix  of  the  interference  is 
given  as 

R,  =  E{(W  V  +  n)(W  V  +  n)"}  =  W"c+c+HW  +  ct,;I  (6) 
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where  <7  ~  is  the  noise  power  and  I  is  the  identity  matrix  and 

the  noise  is  assumed  to  be  zero  mean  white  Gaussian  noise 
(WGN).  Invoking  the  assumption  that  the  large  targets  are 
independent  complex  random  variables  with  variance  defined 

I  +l2 

as  C*  ,  the  processor  based  on  the  optimum  detector  can  be 

written  as  (see  [13]  for  a  more  complete  method  of 
derivation) 

c  =  W(WHc+c+"w  +  o-,;I)^y.  (7) 

This  processor  acts  as  a  matched  filter  while  canceling  the 
previously  detected  large  targets. 

IV.  Strategies  in  applying  CLEAN  Algorithms 

Having  multiple  processors  (i.e.  correlator,  CLEAN 
Deconvolver  and  Reformulated  CLEAN  Detector)  allows  the 
development  of  strategies  to  use  all  these  processors  to  obtain 
the  advantages  of  each.  The  first  algorithm  to  try  is  the 
CLEAN  Deconvolver.  If  the  SNR  is  high  enough,  then  the 
CLEAN  Deconvolver  will  provide  the  best  detection  and 
amplitude  estimation  of  the  multiple  targets  in  the  radar’s 
range  gate.  If  the  smallest  targets  of  interest  have  a  low  SNR 
such  that  they  can  not  be  seen  with  the  CLEAN  Deconvolver, 
then  the  larger  targets  detected  by  the  CLEAN  Deconvolver 
can  be  applied  to  the  Reformulated  CLEAN  Detector  to  allow 
it  to  detect  the  smallest  targets.  If  the  largest  targets  have  such 
a  small  SNR  that  they  can  not  be  detected  by  the  CLEAN 
Deconvolver,  then  they  can  be  detected  with  the  correlator. 
These  targets  can  then  be  applied  to  the  Reformulated  CLEAN 
Detector,  allowing  the  detection  of  the  smaller  targets.  These 
strategies  are  illustrated  in  the  next  section. 


V.  Example  Performance  Results 

The  example  problem  considered  is  that  of  detecting  small 
targets  in  the  presence  of  a  single  large  target.  For  this 
simulation  the  target  Doppler  was  ignored  (set  to  zero).  This 
was  done  for  simplicity  of  presentation.  The  techniques  of 
[10]  have  been  found  to  be  adequate  for  addressing  the 
Doppler  issue  for  this  problem. 

Challenges  addressed  in  this  simulation  include  factional 
spacing  of  targets,  imperfect  knowledge  of  the  PSF  (i.e., 
corrupted  by  noise),  low  SNR  coupled  with  a  high  time- 
bandwidth  pulse.  These  issues  taken  together  severely  limit 
the  performance  of  the  CLEAN  Deconvolver.  The  SNRs 
quoted  are  taken  prior  to  pulse  compression  (ie.,  SNR  of  a 
single  sample). 

The  problem  to  be  solved  in  this  example  is  illustrated  in 
Figure  3.  Figure  3  shows  the  correlator  output  of  a 
hypothetical  radar  using  a  32  chip  derivative  phase  shift 
keying  waveform  with  each  chip  repeated  16  times  (i.e.,  512 
chip  waveform).  The  data  y(t)  is  sampled  twice  per  chip  to 
achieve  the  Nyquist  rate.  Thus,  the  transmitted  waveform  is 
1024  samples  and  the  response  of  the  correlator  is  2047 
samples.  Note  that  for  the  waveform  used  the  range  time 
sidelobes  are  only  17  dB  down  from  the  peak  of  the  target 


response.  There  is  a  large  target  of  35  dB  SNR  in  the  center 
and  three  smaller  ones  with  a  SNR  of  26  dB.  Only  the  large 
target  can  be  reliably  detected  due  to  the  range-time  sidelobes 
of  the  waveform  used.  The  problem  to  be  solved  is  to  detect 
the  three  small  targets  in  the  presence  of  the  large  target. 


Correlator  Output  All  Targets 


detected  dueto  sidelobes. 

Figure  4  shows  the  output  of  the  CLEAN  Deconvolver 
processing  the  same  data.  Here  all  four  targets  can  be  seen. 
Figure  4  shows  the  strength  of  the  CLEAN  Deconvolver, 
namely  that  targets  can  be  localized  to  a  single  range  sample. 
It  should  also  be  noted  that  the  targets  are  not  injected  at  the 
range  cell  boundaries,  but  are  put  in  straddling  range  cells  by 
oversampling  the  input  data  and  then  filtering  and 
downsampling.  Figure  4  then  shows  the  location  of  all  four 
targets  to  nearest  range  sample. 

CLEAN  Deconvolver  Output 

0 - 1 - 1 - - 1 - 1 - 


-10 - 1 - -1 - L 


-70 - 1 - 1 - 1 - 1 - 

0  500  1000  1500  2000  2500 

Range  Cell 

Figure4.  All  FourTargets  are  seen  by  the  CLEAN  Deconvolver. 


The  next  situation  to  be  considered  is  the  case  where  the 
smaller  targets  have  such  a  low  SNR  (12  dB)  that  they  can  not 
be  seen  with  the  CLEAN  Deconvolver.  This  is  illustrated  in 
Figure  5  where  the  noise  and  errors  in  the  estimate  of  the  PSF 
prevent  the  detection  of  the  smaller  targets.  While  not  shown, 
the  correlator  will  also  fail  to  detect  the  smaller  targets. 
However,  for  the  correlator  the  reason  for  non-detection  is  the 
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range -time  sidelobes  of  the  waveform  since  the  SNR  of  the 
small  targets  is  large  enough  to  be  otherwise  detected. 


CLEAN  Deconvolver  Output 


Figure  5.  Only  the  largest  target  is  seen  by  the  CLEAN  Deconvolver. 


Using  the  CLEAN  Deconvolver  to  estimate  the  previously 
detected  target  and  putting  it  into  c+  enables  the  Reformulated 
CLEAN  Detector  to  be  used  to  detect  the  smaller  targets.  The 
Reformulated  CLEAN  Detector  has  the  effect  of  eliminating 
the  large  target  that  was  detected  by  the  CLEAN  Deconvolver. 
This  is  illustrated  in  Figure  6.  Here  the  three  smaller  targets 
are  detected  with  virtually  no  interference  from  the  large 
target.  There  is  however,  some  minor  raising  of  the  time- 
sidelobes  caused  by  the  Reformulated  CLEAN  Detector. 


CLEAN  Detector  Output 


Figure  6.  The  Reformulated  CLEAN  Detector  sees  the  three  small 
targets  while  suppressing  the  large  target. 

The  next  situation  to  be  analyzed  is  the  case  where  the 
larger  target  has  a  SNR  too  low  to  be  detectable  by  the 
CLEAN  Deconvolver.  (The  next  section  quantifies  the 
differences  in  output  SNR  of  these  processors.)  For  this 
example,  the  larger  target  SNR  is  8  dB  while  the  three  smaller 
targets  have  a  SNR  of  -9  dB.  In  Figure  7  we  see  the  output  of 


the  CLEAN  Deconvolver  where  the  large  target  is 
undetectable  because  the  SNR  is  too  low. 


CLEAN  Deconvolver  Output 


-30 


-35 


-45 - 1 - 1 - 1 - 1 - 

0  500  1000  1500  2000  2500 

Range  Cell 

Figure7.  Output  of  the  CLEAN  Deconvolver.  The  SNR  of  the  large 
target  is  too  small  to  allow  detection  by  the  CLEAN  Deconvolver. 


For  the  waveform  used  here  (time-bandwidth  product  of 
512)  it  should  increase  the  SNR  by  27  dB.  Hence  the 
correlator  can  see  the  large  target.  This  is  shown  in  Figure  8 
as  well  as  the  fact  that  the  time-sidelobes  from  the  correlator 
prevent  the  detection  of  the  three  smaller  targets. 


Correlator  Output  All  Targets 


detected. 


Because  the  larger  target  can  be  detected  by  the  correlator 
(see  Figure  8),  the  amplitude  and  position  of  the  target  can  be 
estimated  by  the  correlator.  This  data  is  used  to  populate 
c+  indicating  that  the  large  target  is  now  an  interfering  point 
scatter.  Using  this  information,  the  Reformulated  CLEAN 
Detector  can  detect  the  three  smaller  targets.  This  is 
demonstrated  in  Figure  9. 
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Based  on  this  the  SNR  for  the  Reformulated  CLEAN 
Detector  for  the  Vth  range  cell  is 


CLEAN  Detector  Output 


targets  can  be  detected  once  the  large  target  is  removed. 

While  the  examples  in  this  section  only  had  one  large  target, 
these  techniques  could  be  used  for  dealing  with  any  number 
of  large  interfering  targets,  as  long  as  they  can  be  individually 
detected  and  removed.  In  fact  these  techniques  work  even 
better  as  the  targets  are  larger  having  a  large  SNR. 

VI.  SNR  Comparisons 

Since  is  it  is  well  known  that  the  matched  fdter  produces 
the  maximum  output  SNR  for  the  situation  of  WGN  this  begs 
the  question  as  to  what  is  the  SNR  performance  of  the 
CLEAN  Deconvolver  and  the  Reformulated  CLEAN 
Detector.  In  this  section  the  SNR  for  CLEAN  processors  will 
be  determined. 

First  for  the  CLEAN  Deconvolver  the  signal  output  due  to  c 
can  be  determined  from  (1)  and  (2)  as 

c  =  (WW"r‘WW"c  =c.  (8) 

Thus  the  signal  output  power  is  |c|2.  The  output  noise 
correlation  matrix  using  (2)  and  noting  that  n  is  zero  mean 
white  Gaussian  noise  is 

R„  =  E{(WW“ )~‘  WnnH  Wa  (WWfl )  1 }  =  (t,!(WW“  )  1 .  (9) 

This  gives  the  SNR  for  the  Cth  range  cell  as 

SNRc_Decon=|c,|2/a„2(WW  Hy(lk).  (10) 

Next  the  Reformulated  CLEAN  Detector  output  due  to  the 
signal  is  determined  from  (1)  and  (7)  as 

c  =  W(W''cVi,W  +  (T,;irlW,'c.  (11) 

The  output  noise  correlation  matrix  for  the  Reformulated 
CLEAN  Detector  is  determined  as 

R„  =E{W(W"c+c+HW  +  CT,(lr1nn"(WHc+c+"w  +  CT,;lr1WH}(12) 
This  can  be  simplified  as 

R  =  CT,fW(W"cV"w  +  CT,;IL1(W"cV"w  +  <7,;irlWfl.  (13) 


SNR  c_Det  \ck  |  /  (14) 

where  c  and  R„  are  defined  in  (11)  and  (13)  respectively. 

The  ratio  of  the  SNR's  for  the  two  CLEAN  processors  to 
the  matched  filter  is  given  in  Figure  10.  Figure  10  gives  the 
relative  performance  for  the  example  calculations  in  the 
previous  section  with  one  large  target.  Thus  it  can  be  seen  that 
the  Reformulated  CLEAN  Detector  suffers  almost  no  loss 
except  in  the  vicinity  of  the  large  target  that  it  is  suppressing. 
In  contrast  the  CLEAN  Deconvolver  suffers  significant  SNR 
loss  compared  to  the  matched  filter.  It  should  be  noted  that 
the  loss  is  waveform  dependent  and  is  likely  to  be  large  for  the 
large  time -bandwidth  product  (i.e.,  512)  waveforms  used  in 
the  previous  section. 


SNR  Ratios 


Figure  10.  This  figure  shows  the  SNR  disadvantage  the  CLEAN 
Deconvolver  and  the  CLEAN  Detector  have  compared  with  the  Correlator 


VII.  Comparisons  to  Other  Approaches 

At  the  suggestion  of  the  reviewers  this  section  compares 
and  contrasts  the  approach  in  this  paper  to  two  others.  The 
first  is  [8]  which  applies  CLEAN  to  Space  Time  Adaptive 
Processing  (STAP)  problem.  The  authors  apply  the  two 
dimensional  CLEAN  algorithm  (much  like  radio  astronomers 
[2])  to  resolve  targets  and  clutter.  The  authors  use  Taylor 
weighting  to  lower  the  sidelobes  of  the  PSF  to  aid  in  the 
function  of  the  CLEAN  algorithm  to  resolve  mainlobe  and 
sidelobe  responses.  The  use  of  this  data  is  to  remove  array 
calibration  errors.  Thus,  the  measure  of  performance  is  not  in 
detecting  all  the  clutter  or  target  sources  but  in  how  well  the 
output  of  CLEAN  allows  the  data  cube  to  be  calibrated  such 
that  subsequent  processing  can  reliably  detect  targets. 
Probably  the  biggest  difference  in  this  approach  and  this  paper 
(besides  that  [8]  is  a  two  dimensional  problem)  is  that  the 
sidelobes  of  the  PSF  are  controlled  to  limit  sidelobes.  The 
PSFs  in  this  paper  had  range-time  sidelobes  on  the  order  of  17 
dB  and  one  of  the  strengths  of  this  approach  is  to  deal  with 
PSFs  with  high  range-time  sidelobes. 
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The  next  paper  to  compare  is  [11]  which  uses  a  technique 
called  Adaptive  Pulse  Compression  (APC).  The  APC 
approach  it  to  modify  the  correlator  based  on  nearby  targets  to 
minimize  the  mean  square  error  in  estimating  the  range  cell  of 
interest.  The  statistical  approach  of  [11]  is  Bayesian  while 
that  of  this  paper  is  a  classical  approach.  Thus  [11]  assumes 
that  c  is  a  random  vector  of  known  distribution  and  seeks  to 
minimize  the  mean  error  in  estimating  it.  In  contrast  the 
CLEAN  Deconvolver  assumes  that  c  is  deterministic  and 
produces  an  unbiased  estimate  of  c  that  has  an  error  variance 
that  is  lower  than  all  other  estimators.  The  CLEAN 
Deconvolver  does  not  have  to  estimate  any  intermediate 
values  to  produce  its  estimate  of  c  as  is  the  case  of  APC.  The 
Reformulated  CLEAN  Detector  on  the  over  hand  is  derived 
using  detection  theory  as  opposed  to  estimation  theory.  To 
apply  it  one  must  independently  estimate  c+  (vector  of  large 
targets).  The  Reformulated  CLEAN  Detector  can  also  be 
described  as  a  whitening  matched  filter.  In  addition,  its 
optimality  claim  is  that  it  produces  a  test  statistic  that 
maximizes  the  probability  of  detection  for  a  given  false  alarm 
probability  as  long  as  n  and  c+  are  jointly  Gaussian  distributed. 


VIII.  Summary  and  Conclusion 

This  paper  extended  the  previously  developed  CLEAN 
Deconvolver  and  CLEAN  Detector  by  developing  the 
Reformulated  CLEAN  Detector.  Based  on  these  processors 
and  the  correlator,  strategies  were  described  as  to  how  these 
processors  could  be  combined  to  address  the  problem  of  low 
SNR  targets.  The  potential  performance  of  this  approach  was 
then  given  in  section  V  where  it  was  shown  that  small  targets 
could  be  detected  in  the  presence  of  large  targets. 
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